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ABSTRACT 

We consider the partition function for a matrix model with a global 
unitary invariant energy function. We show that the averages over the par- 
tition function of global unitary invariant trace polynomials of the matrix 
variables are the same when calculated with any choice of a global unitary 
fixing, while averages of such polynomials without a trace define matrix- 
valued correlation functions, that depend on the choice of unitary fixing. 
The unitary fixing is formulated within the standard Faddeev-Popov frame- 
work, in which the squared Vandermonde determinant emerges as a factor 
of the complete Faddeev-Popov determinant. We give the ghost representa- 
tion for the FP determinant, and the corresponding BRST invariance of the 
unitary-fixed partition function. The formalism is relevant for deriving Ward 
identities obeyed by matrix-valued correlation functions. 
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Over the years there has been considerable interest in matrix models from various 
points of view. Matrix models are used to approximate quantum many body systems and 
quantum field theories [1], and have deep connections with string theories [2]. They also 
have been studied as classical statistical mechanical systems, from which quantum behavior 
emerges under certain conditions [3]. A common issue that arises in all of these applications 
is dealing with an overall global unitary invariance transformation of the partition function. 
Typically, in matrix model calculations this overall invariance is partially integrated out as 
a first step, thus eliminating a U(N)/U(l) N subgroup of the global unitary group. Our aim 
in this paper is to proceed in an alternative fashion, by using the Faddeev-Popov framework 
to impose a set of unitary invariance fixing conditions, that completely break the SU(N) 
subgroup of the global unitary invariance group U(N). One can think of our construction 
as a type of polar decomposition, based on modding out the action of the SU(N) subgroup. 
This allows one to define matrix- valued correlation functions, which give additional structural 
information about the system, but which (like gauge potentials in gauge field theory) depend 
on the choice of unitary fixing. A complete global unitary fixing is needed for the application 
of matrix models to emergent quantum theory developed in Ref. [3], so as to be able to 
construct matrix ensembles that do not integrate over the spacetime translation group of 
the emergent theory. The formalism that we develop here may well find other matrix model 
applications as well. 

Let Mi, ...,M D be a set of N x N complex self-adjoint matrices, and let us take as 
the energy function 

U[{M}]=TrH(M 1 ,...,M D ) , (1) 

with H a self-adjoint polynomial in its arguments constructed using only c-number coeffi- 
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cients (i.e., no fixed, non-dynamical matrices appear as coefficients in constructing H). Then 
the corresponding partition function Z is defined by 

dM exp(-H) , (2) 

with 

D 

dM = ]Jd[M d } , (3a) 

d=l 

and with the integration measure d[M] for the self-adjoint matrix M defined in terms of the 
real and imaginary parts of the matrix elements My of M by 

d[M] = Y[dMa YldReMijdLmMij . (36) 

i i<j 

As is well known, the measure d[M] is unitary invariant, in other words, if U is a fixed N x N 
unitary matrix, then 

d[U*MU] = d[M] . (4a) 

If we make the same unitary transformation U on all of the matrices Md , d — 1, D, then 
by our assumption that H involves no fixed matrix coefficients, H is invariant by virtue of 
the cyclic property of the trace, 

H[{U j MU}] = H[{M}] . (46) 

Thus, Eqs. (4a) and (4b) together imply that the partition function Z has a global unitary 
invariance. 

The global unitary invariance of Z must be taken into account in calculating corre- 
lations of the various matrices M d averaged over the partition function. Let Q[{M}] be an 



arbitrary polynomial in the matrices Mi, Mp constructed using only c- number coefficients, 
so that under global unitary transformations, Q transforms as 

Q[{U ] MU}\ = U ] Q[{M}\U . (5a) 

Correspondingly, let 

Q = TrQ , (56) 

so that Q is a global unitary invariant. One can now consider the calculation of averages 
of Q and of Q respectively over the ensemble defined by Eq. (2). In the case of the trace 
polynomial Q one has 

(Q)av = Z- 1 J dM exp(-H)Q , (6a) 

which because of the global unitary invariance involves an overall structure-independent 
unitary integration that is typically done as the first step, by using Mehta's change of vari- 
ables [4] for one of the matrix arguments on which Q depends. Let us now consider the 
corresponding average of the polynomial Q over the ensemble, 

(Q)av = Z- 1 J dMexp(-H)Q . (66) 

Making a global unitary transformation on all of the matrix integration variables, and using 
the invariance of dM and of H given in Eqs. (4a,b), and the covariance of Q given in Eq. (5a), 
we then find that 

(Q)av = UHQ) ay U , (7a) 

for all unitary matrices U. Thus by Schur's lemma (which applies since U (N) acts irreducibly 
on the complex N dimensional vector space) (Q)av must be a c-number multiple of the unit 



matrix, so that by taking the trace, we learn that 



(Q) AV = iV- 1 (Q) AV 



(76) 



and all nontrivial matrix information (e.g., the unitary orientation and nontrivial operator 
properties) contained in Q has been lost. 

In order to retain access to the matrix information contained in Q, let us then proceed 



global unitary invariance group has been fixed. (The full global unitary invariance group is 
the product of this SU (N) with a global U(l) that is an overall phase times the unit matrix; 
since this U(l) commutes with Q, averaging over it causes no loss of the matrix information 
contained in Q, and so fixing the overall U(l) is not necessary.) We then define the average 
of Q over the unitary-fixed ensemble as 



the partition function in which the global unitary invariance has been broken, and an orienta- 
tion on the iV-dimensional vector space has been fixed. Clearly, the procedure just described 
is a global unitary analog of the gauge fixing customarily employed in the case of local gauge 
invariances. If we change the recipe for fixing the global unitary invariance, then the average 
defined by Eq. (8a) will change in a manner that is in general complicated. However, we will 
show that the average of Q in the unitary-fixed ensemble is independent of the fixing and is 
equal to that defined in Eq. (6a) by averaging over the original ensemble, so that 



in an alternative fashion. Let us define a measure dM in which the SU (N) subgroup of the 




(8a) 



with 




(86) 



(Q)av = (Q)av 



(9) 
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In other words, the average of the trace of Q takes the same value for any choice of unitary- 
fixing. To make an analogy with local gauge fixing in gauge theories, the trace polynomials Q 
are analogs of gauge invariant functions, while polynomials Q without a trace are analogs of 
gauge-variant quantities. Just as the gauge-variant potentials contain useful information in 
gauge theories, the unitary fixing- variant averages of polynomials Q contain useful structural 
information about matrix models. 

To prove Eq. (9), we proceed by analogy with the standard Faddeev- Popov procedure 
used for local gauge fixing. Let us write an infinitesimal SU (N) transformation in generator 
form as U — exp(C), with G anti-self- adjoint and traceless. We take as the N 2 — 1 infinites- 
imal parameters of the SU(N) transformation the real numbers gj , j — 1, TV 2 — 1, 
with those for j = 1,...,N(N — 1) given by the real and imaginary parts of the off- 
diagonal matrix elements of G, that is, by ReG^ and ImGV, for i < j. The remaining 
ones for j = N(N — 1) + 1, ...,N 2 — 1 are given by the differences of the imaginary parts 
of the diagonal matrix elements of G, that is, by Im(Gn — G22 ),•••, Im(Gn — Gnn)- Let 
fj({M}) , j = 1, N 2 — 1 be a set of functions of the matrices M 1 , M D with the prop- 
erty that the equations fj({M}) = , j = 1,...,N 2 — 1 completely break the SU(N) 
invariance group, so that the only solution of fj({M + [G, M]}) — , j — 1, ...,N 2 — 1 is 
9 j — , j = 1, N 2 — 1. We consider now the integral 



Here Q is a global unitary invariant function of the matrices Mi, Mo, such as a trace poly- 
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(10a) 



with the function taken as 



N 2 -l 



Km] = n 



(106) 



nomial Q or any function of trace polynomials (for example the partition function weight 
exp(— H)). Equation (10a) has the standard form of the Faddeev-Popov analysis, as for- 
mulated for example in the text of Weinberg [5] (except that when one is dealing with a 
non-compact local gauge invariance, where the limits of integration lie at infinity, one can 
take the function K to be a general function of gauge variant functions /_,•; in the compact 
case considered here, the delta functions of Eq. (10b) must be used in order to make the 
integration limits irrelevant.) The standard FP argument then shows that the integral in 
Eq. (10a) is independent of the constraints fj. Briefly, the argument proceeds by replac- 
ing the dummy variable of integration dM by dM v , where M v = V'MV, and integrating 
over the SU (N) matrix V. The group property of unitary transformations together with 
the chain rule then converts the determinant in Eq. (10a) into a Jacobian transforming the 
V integration into an integration over the constraints fj, permitting the delta functions in 
Eq. (10b) to be integrated to give unity. This shows that the result is independent of the 
constraints, and that it is the same as the result obtained by integrating over the original 
unfixed ensemble, thus establishing Eq. (9). Clearly, this argument works only when the 
function Q is a unitary invariant, so that it has no dependence on V. For example, if Q 
is replaced by a polynomial in the matrices without an overall trace, then the unitary fix- 
ing constraints cannot be eliminated by integrating over V, and the result depends on the 
unitary fixing in a complicated way. 

A specific realization of the general unitary fixing can be given when D > 2, so that 
the set of matrices Mi,..., M D contains at least two independent self-adjoint matrices A — Mi 
and B = M 2 . We take the functions fj , j = 1, N 2 — 1 to be linear functions of A and 
B, constructed as follows. As the fj for j = 1, N(N — 1) we take the real and imaginary 
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parts of the off-diagonal matrix elements of A, that is, the functions Re A^ and Im A^ for 
% < j. Equating these functions to zero forces the matrix A to be diagonal. The N remaining 
diagonal unitary transformations then commute with A, so that no further conditions can be 
furnished by use of A alone. However, the diagonal SU(N) transformations can always be 
used to make the off-diagonal matrix elements in the first row of the second matrix B have 
vanishing imaginary parts, leaving a residual Z^~ l symmetry that is broken by requiring 
these matrix elements to have positive semidefinite real parts. So for the remaining conditions 
j) for j = N(N — 1) + 1, N 2 — 1, we take the N — 1 functions ImB^ , j > 1, and we 
restrict the integrations over KeBij , j > 1 to run from to oo. Since the function K chosen 
in Eq. (10b) enforces the conditions fj = in a sharp manner, they can be used to simplify 
the expression for the Faddeev-Popov determinant. A simple calculation now shows that 
when the fj all vanish, the matrix elements of the commutator [G, M] needed in Eq. (10a) 
are given by 



with R a remainder containing only off-diagonal elements G^j of the matrix G. Since 
Eq. (11a) shows that the matrix 



Re[G, A]ij =ReG ij (A jj - An) 



lm[G, A]ij =\m.G ij (A jj - A u ) 



(11a) 



Im[G, B] X j =ReS y Im(Gn - Gjj) + R 




(116) 



is triangular (its upper off-diagonal matrix elements are all zero because R has no dependence 
on the diagonal matrix elements of G), its determinant is given by the product of its diagonal 
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matrix elements. Thus we have 



A S de t ( a/ - ({M+|G - Ml » i ) =n<*.-^)'nitafl« • < i2a > 

V • ' ' i<j 3=2 

the first factor of which is the familiar squared Vandermonde determinant. Substituting 
Eqs. (10b) and (12a) into Eq. (10a), we thus arrive at the formula for the unitary-fixed 
integral 

/D N N 

HdiMdliYldAudB^illdReB^i J[ tflteB o -tflm£ o -)A0[{M}] , (126) 

d=3 i=l j=2 2<i<j 

with the integrals over ReBij , j = 2,...,N in Eq. (12b) running over positive values only. 
The part of this analysis involving only a single matrix A is well-known in the literature [6]; 
what has been added here is the complete SU (N) fixing obtained by imposing a condition on 
a second matrix B as well. The part of Eqs. (12a,b) involving each B±j is just a planar radial 
integral / °° pdp, with p = \Bij\ = ReBij, where the associated angular integral J Q 27r d<p has 
been omitted because it corresponds to a U(l) factor that has been fixed by the condition 

= 0. 

With this choice of unitary fixing, the unitary fixed average Q = (Q) A y defined in 
Eq. (8a) has a characteristic form that is dictated by the symmetries of the unitary-fixed 
ensemble. Since the unitary fixing conditions are symmetric under permutation of the basis 
states with labels 2,3, ...,N, and since this permutation is also a symmetry of the unfixed 
measure dM, the matrix Q must be symmetric under this permutation of basis states. Thus, 
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there are only five independent matrix elements, 



Q 



a 



■P, j 2 V 



7 , J 



N 



(12c) 



Qn 



5 , i = 2,..., 



iV 




e,2<i^j<N 



In this notation, the original unfixed average Q = {Q)av defined by Eq. (7b) is given by 



showing explicitly that there is a loss of structural information in using the unfixed av- 
erage. But even the unitary-fixed average has a structure that is greatly restricted as 
compared with a general N x N matrix. (Similar reasoning applies to the partial uni- 
tary fixing in which one only imposes the condition that A should be diagonal. Since this 
condition is symmetric under permutation of the basis states with labels 1, N, the par- 
tially unitary fixed average of a polynomial Q defined by integrating with the measure 
(IL>=i dAii) ( Ilj<j(Ai — Ajj) 2 ) Yld=2 d[Md] must also have this permutation symmetry, and 
thus must be a c-number times the unit matrix.) 

We now introduce ghost integrals to represent the determinant A. Let u>ij and cuij 
be the matrix elements of independent N x N complex anti-self-adjoint Grassmann matrices 
uj and Cj. We take uj to be traceless, Trcj = 0, while we take Cj to have a vanishing 11 matrix 
element, Wu = 0. The integration measure for ui is defined by 




(12d) 



N 

du = Y\ dRecUijdlmcUij Y\ d\m.{ujj — Wn) 

i<j j=2 



(13a) 
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while the integration measure for u is taken as 

N 

du — Y\ dReuijdlmuij J^J dlmujj . (136) 

i<j j=2 

We can now use these Grassmann matrices to give a ghost representation of the factors in 
Eq. (12a) involving the matrices A and B. Since the matrix A is diagonal, we have 

Trcifcj, A] = y^^ujjijAu - Ajj)uij . (14a) 

Hence up to an overall sign, the square of the Vandermonde determinant rii<j(Ai — ^-jj) 2 
is given by the ghost integral 

J d'ud'Cu exp(Tiu[u, A]) , (146) 

with the diagonal factors dlm(ujj — Wn), dlmujjj , j — 2, ...,N omitted from the primed 
integration measures d'uj and d'Cij. To represent the second factor in Eq. (12a) as a ghost 
integral, we use the diagonal matrix elements of u and Co in an analogous fashion. Thus, up 
to a phase, the factor n^=2 R- e -Sii is given by the ghost integral 

/N N 
dlm(ujjj — u>ii)dlmCujj exp ( Cujj(ReBij)i(ujj — Wn)) . (14c) 

3=2 3=2 

By defining a matrix X by Xu = 0; Xij = , 2 < i,j < N; Xij = Xji = ^Cjjj(ujjj — 
Uu) , j = 2, ...,iV, the exponent in Eq. (14c) can be written as TrX_B, so that Eq. (14c) 
becomes 

N 

/ IT — UiijdlmUjj exp(TrX_B) . (14o?) 

Combining Eqs. (13a,b) and (14b,c,d), we see that up to an overall phase the determinant 
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A introduced in Eq. (12a) has the equivalent ghost representations 



N 



A oc / dudw exp (Trcufcu, A] + ^^ujjj(ReB 1 j)i(u!jj — oon)) 

J i=2 (15) 

oc J dcudujexp(Tru;[uj,A\+TrXB) . 
Yet another equivalent form is obtained by noting that 

[B.u^j = Bijlujj - u n ) + S , (16a) 

with the remainder S denoting terms that only involve matrix elements uiij with % ^ j. The 
remainder S makes a vanishing contribution to the Grassmann integrals when Eq. (16a) is 
substituted for B\ 3 i{uj 3] - — wu) in Eq. (15), since one factor of (cujj -wu) for each j = 2, N 
is needed to give a nonvanishing integral, and each such term in the exponent is already 
accompanied by a factor Ujj, so that terms with additional such factors vanish inside the 
Grassmann integrals. (We are just using here the fact that with (, ( Grassmann variables, 
J d(d(exp((W( + U() = W, with no dependence on U .) Since the diagonal matrix elements 
of uj are pure imaginary, Eq. (16a) implies that 

{ReBij)i{ujjj — cuu) = —lm[B,cu]ij + ImS* , (166) 

which when substituted into Eq. (15) gives the alternative formula 

r N 

A oc / dudu exp (Ttu[u, A] —^^Ujjlm[B,uj]ij) . (17) 

This formula will be used to establish a BRST [7] symmetry, the topic to which we now turn. 

To formulate a BRST invariance transformation corresponding to Eq. (17), we rewrite 
the product of 5 functions in Eq. (10b) and the half-line restriction on the integrals over KeBij 
in terms of their Fourier representations, by introducing three sets of Nakanishi-Lautrup [8] 
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variables. One set are the elements hij of a self-adjoint N x N matrix h with vanishing 
diagonal matrix elements, so that ha — , i — 1, N. The integration measure for this set 
is defined as 

dh — JT dRehijdlmhij . (18a) 

i<j 

The second set are N — 1 real numbers Hj , j — 2, N, with integration measure 

N 

dH = Y[dH j . (186) 
i=2 

In terms of these variables, the product of 5 functions of Eq. (10b) can be represented (up 
to an overall constant factor) as 

N 2 -l . N 

Y[ 6{fj) oc / dhdH exp{iTihA + ij^ H M B ij) ■ ( 19a ) 

3=1 J 3=2 

The third set are — 1 complex numbers kj , j = 2,...,N, integrated along a contour on 
the real axis with integration measure 

N 

dk = Y[dk j /(k j -ie) , (196) 

i=2 

with infinitesimal positive e. These can be used to insert a product of step functions 

Uf=2 e ( ReB i3) into E q- ( 12b )> 

N „ N 

dkexp (i kjReBij) , (19c) 

3=2 J 3=2 

allowing the integrals over the ReSy in Eq. (12b) to be taken from — oo to oo. 

Defining a matrix Y by Y u = 0; Y i3 r = , 2 < i,j < N; Y Xj = -Y n = -\Hj, the 
second term in the exponent in Eq. (19a) can be rewritten as i Yl!j=2 HjlraBij = TrYB, and 
so an alternative form of Eq. (19a) is 

N 2 -l 

(19d) 



Y[ S(fj)(x f dhdHexp(iTrhA + TrYB) 

3=1 J 
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Similarly, defining a matrix Z by Z\\ = 0; Z i3 - = , 2 < i,j < N; Z\ 3 - = Z 3 \ = \ik 3 , the 
exponent in Eq. (19c) can be rewritten as J2f=2 ik 3 ReB l3 = TrZB, and so an alternative 
form of Eq. (19c) is 



N 

Jj0(ReSij) oc / dkexp(TrZB) . (19e) 

j=2 J 



These equations allow us to write Eq. (12b) in terms of the unrestricted measure dM, and 
the ghost representation of A, as 



J =C I dMdhdHdkdudCj 

N 

I 

i=2 



/■ 

x exp (iTihA + Ttlj[u, A] + ^{iHAm.B l3 + ikjReB l3 - d)^Im[S, uo] l:j ))g[{M}] 

i=2 

=C J dMdhdHdkduduexp (iTrhA + Tiu[uj, A] + Tr(X + Y + Z)B)g\{M}} , 

(20) 

with C an overall constant factor. The first representation of J in Eq. (20) will be used to 
establish a BRST invariance, while the second will be used to discuss Ward identities obeyed 
by the matrix-valued correlations. 

We now show that the first integral in Eq. (20) is invariant under the nilpotent BRST 



15 



transformation 

5A=[A,u}6 , 
5B=[B,uj}6 , 
5M d =[M d ,u}9 , d = 3,...,D , 
5uj =uj 2 9 , 

SuJij = - ih i:j e , % ^ j , (21) 
Sujjj = — iHjO ,j = 2, iV , 

5h =0 , 
5Hj =0 , 
5kj =0 , 

with 8 a c-number Grassmann parameter. (The part of this transformation involving uj is 
patterned after the BRST transformation for the local operator gauge invariant case studied 
by Adler [9].) We first remark that since Eq. (21) has the form of an infinitesimal unitary 
transformation with generator ujQ acting on the matrix variables M d , the global unitary 
invariant function C?[{M}] and the matrix integration measure dM are both invariant. We 
consider next the terms in the exponent in Eq. (20). From Eq. (21) we have 
5[A, uj\ =[5 A, u] + [A, 5u] = [[A, u]0, uj] + [A, u 2 6] 

(22a) 

= - (w[A, uj] + [A, uj]uj)9 + [A, lu 2 ]9 = -[A, uo 2 ]6 + [A, uj 2 ]9 = . 
Hence for the terms in the exponent of Eq. (20) involving A, we get (using the fact that A 
is diagonal) 

5(tTrhA + TtCu[uj, A]) = iTrhSA + Tt(5Co)[uj, A] 

(226) 

= iTrh[A,u]e + Tr(-ihe)[u,A]=0 . 
For the terms in the exponent of Eq. (20) involving B but not involving the parameters kj, 
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inside the summation over j we have 



5(iHjlmBij — Ujjlm[B,u]ij) = iHjlmSBij — (5ujjj)lm[B,uj]ij 



(22c) 



iHj1m[B, u^xjO + iHj8Im[B, w] y = 



since 5[B,uj] = by the same argument as in Eq. (22a). So the entire exponent of the first 
representation in Eq. (20) is BRST invariant, apart from the kjReBij terms. But the shifts 
in the KeBij are linear in uj while not involving Co. Thus (since we shall see shortly that 
the integration measures are invariant), the shifts in the terms in the exponent involving the 
products kjReBij make a vanishing contribution to the Grassmann integrals, by an argument 
similar to that used to justify the neglect of S in Eq. (16a). 

An alternative method of including the step functions, that leads to a manifestly 
BRST invariant integrand, is to include in the exponent in the first representation of Eq. (20) 
an additional term — ^2f =2 KjKe[B, with Kj auxiliary Grassmann parameters that are 
not integrated over. This term is linear in uj but does not involve u, and so again makes a 
vanishing contribution to the Grassmann integrals in Eq. (20). The BRST transformation 
of Eq. (21) is then augmented by the rule 5kj = —ikj9, with the result that the combination 
ikjReBij — KjRe[B, u>]ij is manifestly BRST invariant. 

Continuing the BRST analysis, since Tiar = — Trrcr for any two Grassmann odd 
grade matrices r and a, we have Tru; 2 = — Tru; 2 = 0, and so the condition that uj should 
be traceless is preserved by Eq. (21). (On the other hand, u)\ x is nonzero even when un is 
zero, which is why we must use a traceless condition, rather than a condition Wn = 0, for 
uj.) Also, letting * denote complex conjugation, since 



the property that uj is anti-self-adjoint is preserved by Eq. (21). The integration measures 
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dh and dH are trivially invariant, while the measure du is invariant because 5u has no 
dependence on Cj. Since 

S(dcOij) = d(5cu)ij = d(u 2 9)ij = (ujduj + (dcu)uj).£ , (23a) 

we have 

S(duij) = (uJudujij + duJijUJjj)9 + .... = duij(ujj - uou)9 + (236) 

with ... denoting terms that contain only matrix elements du^ y with ^ Hence 

there is no Jacobian contribution from the diagonal terms in the measure du, while the 
Jacobian arising from transformation of the off-diagonal terms in du> differs from unity by a 
term proportional to 

y^jj - uJii)9 = , (23c) 

and so the measure du is also invariant. Finally, nilpotence of the BRST transformation 
follows from Eq. (22a), and its analogs with A replaced by B or by a general M d) together 
with 

5lo 2 = {5lo, lo} = {u 2 9, u} = u 2 {9, lo} = . (23d) 

This completes the demonstration of the BRST transformation for the first representation 
in Eq. (20). 

The second representation in Eq. (20) can be used to derive Ward identities from 
unitary- fixed expectations of trace polynomials Q; these Ward identities play a central role 
in the arguments for an emergent quantum theory given in Ref. [3]. Employing the specific 
unitary fixing of Eq. (20) in the definition of Eqs. (8a,b), as applied to Q = TrQ, and using 
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the cyclic property of the trace to rewrite Ttlu[lu, A] as Tr{u), u}A, we have 

Z(Q) AV = j d M dhd H dkd ^ mh + + TV(X + K + Z)B) e*p(-H)Q , 

(24) 

with Z given by the expression on the right hand side of Eq. (24) with Q replaced by unity. 
Ward identities follow from the fact that the unrestricted measure dM is invariant under a 
shift of any matrix M d by a constant 5M d , which under the assumption that surface terms 
related to the shift vanish, implies 

= J dMdMHdkdudu>6 Md (exp (Tr[(ih + {Co, to}) A] + Tr(X + Y + Z)B) exp(-H)Q) . 

(25a) 

When H and Q are varied with respect to Md, the factor SMd can be cyclically permuted 

to the right in each term of the varied trace polynomials, giving the formulas 

SB. 
SM t 
5Q 



(256) 

5 Md Q=Tri^-5M d , 



SM d 

which [3] define the variational derivatives of the trace polynomials with respect to the 
operator M d . Carrying through the variations of all terms of Eq. (25a), and dividing by Z, 
we are left with an expression of the form 

= Ti(W d )^5M d . (26a) 

However, since 5M d is an arbitrary self-adjoint matrix, the vanishing of the real and imagi- 
nary parts of Eq. (26a) implies the matrix identity 

= (^) AV • (266) 
19 



For d — 3, D, the variation 5M d in Eq. (25a) acts only on the product exp(— H)Q, and we 
have 

However, for d — 1 and d = 2, corresponding to M 1 = A and M 2 = B, there are additional 
contributions to the Ward identities arising from variations of the traces involving A and B 
in the first exponential on the right hand side of Eq. (25a), which arose from the unitary 
fixing procedure. Explicitly, we have 

Wi=(i/i + {w,w})Q+^-Q^5 , 

(276) 

M + y + z )Q + ^-Q§ . 

Hence from Eq. (20) we are able to get explicit forms of all of the Ward identities, including 
those obtained by varying the matrices singled out in the unitary invariance fixing. Note that 
were we to employ the original ensemble average of Eq. (6a), which has no unitary fixing, 
in deriving the Ward identities, then Eq. (7b) implies that we would only obtain the trace 
of the matrix relation of Eq. (26b). In other words, unitary fixing is essential for extracting 
the full content of the Ward identities; without it, all nontrivial matrix structure is averaged 
out. 
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